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Abstract 

The DBI galileons are a generalization of the galileon terms, which extend the internal galilean 
symmetry to an internal relativistic symmetry, and can also be thought of as generalizations of 
DBI which yield second order field equations. We show that, when considered as local modifica- 
tions to gravity, such as in the Solar system, there exists a region of parameter space in which 
spherically symmetric static solutions exist and are stable. However, these solutions always exhibit 
superluminality, casting doubt on the existence of a standard Lorentz invariant UV completion. 
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I. INTRODUCTION 



The Dirac-Born-Infeld (DBI) action, which describes the dynamics of a brane embedded in 
a higher dimensional spacetime, has provided an important setting within which to study 
inflation [TJ [2], late-time cosmic acceleration [3], tunneling [I], and exotic topological de- 
fects [SHU]. The DBI action has been extensively studied in recent years, and its rather 
special properties are now well-understood. 

At the same time, increasing attention has been paid to induced gravity theories, such 
as the Dvali-Gabadadze-Poratti (DGP) model [TU]. In these, branes in extra dimensions 
again form the basic objects, but the Einstein-Hilbert action for gravity is written both 
in the bulk and on the the branes themselves, leading to a highly nontrivial behavior of 
the resulting 4-dimensional effective theory on the brane. The theory admits a limit which 
contains a scalar field 7r which interacts through a higher derivative cubic coupling possessing 
an internal galilean invariance 5tt = u^x^ + e (with and e infinitesimal constants) and 
second order field equations [TTJ H2] . These general properties can be generalized to higher 
order interactions, and the resulting theories are know as galileons [13] . 

It has recently been shown [H] that the galileon and DBI theories are intimately related. 
Beginning from a co-dimension one probe brane in a 5D Poincare symmetric bulk, a 4D 
action can be formed using solely the 4D Lovelock invariants and the boundary terms of 
5D Lovelock invariants. The resuting action then consists of DBI galileon terms. These 
terms are generalizations of the square root DBI action, in the sense that they share its 
symmetries and yield second order field equations. In a small field limit, the relativistic 
symmetry stemming from the 5D Poincare symmetry and brane reparametrization invariance 
reduces to the galilean symmetry of the DGP model and the DBI galilean terms reduce to 
the galilean terms catalogued in 1 1 3 J . Several authors pT5HT8] have now demonstrated that 
a natural generalization of these galileon models to co-dimension greater than one exists, 
and that many features, including a powerful non-renormalization theorem and a consistent 
effective field theory [16], hold in general. 

In this paper we find general static spherically symmetric solutions of DBI galileon theories, 
and explore their stability. Such an analysis was performed for the ordinary galileons in [13] , 
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and for multi-galileon theories in JT9J 120]. In the case of the DGP model [13], it was found 
that for some choices of parameters, stable solutions exists but always contain superluminal 
signal propagation. We follow the same approach here, extending the results to the DBI 
galileons, and reach similar conclusions. The analysis is only valid in the M p — > oo limit. 
As shown in [21] , the stability of these theories depends on terms suppressed by the square 
of the Planck mass. 



II. DBI GALILEAN TERMS AND EQUATIONS OF MOTION 

We are interested in the generalizations of DBI discussed in [TJ]. The relevant theory consists 
of a single scalar 7i~, in 3+1 dimensions, with an action invariant under the internal relativistic 
symmetry 

5n = - u^-Kd^-K + e , (1) 

with x M the spacetime coordinate, a constant infinitesimal vector, and e an infinitesimal 
constant. 

To construct the action for 7T, one follows the prescription of [TJ] (see also [IB] for further 
details). Consider an embedding of a 3-brane in flat 5-d Minkowski space, X A (x^) (where A 
is the 5-d bulk index and /i the 4-d world- volume index) , and a world- volume action which is 
invariant under world-volume reparametrizations and bulk Poincare transformations. The 
reparametrization invariance forces the action to be a diffeomorphism scalar constructed 
out of the induced metric g^ix) = j^^tGab (X(x)), where Gab is the bulk metric as a 
function of the embedding variables X A (x). Poincare invariance requires the bulk metric to 
be the flat Minkowski metric Gab (A) = tjab- We then fix the gauge X^(x) = and let 
the unfixed degree of freedom be X 5 = 7T, so that the induced metric becomes 

g^u = r),jj, + d^dv-K . (2) 

Any action which is a diffeomorphism scalar, evaluated on this metric, will yield an action 
for 7r having the invariance ([T]), in which is a boost in the fifth direction along with 
a compensating gauge transformation to maintain the gauge choice A M (x) = x M . The 
parameter e, the shift on n, is the translation in the fifth dimension. In addition, the action 
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will have the usual 4-dimensional spacetime Poincare invariance, which combines with the 
boost and the shift ([!]) to form the full 5-d Poincare group. 

The ingredients available to construct such an action are the induced metric g^, the covari- 
ant derivative V /t compatible with the induced metric, the Riemann curvature tensor R p atiu 
corresponding to this derivative, and the extrinsic curvature K^ v of the embedding. Thus, 
the most general action is 



(3) 



For example, the DBI action arises from 

J d 4 x ^~g -> J d 4 x y/1 + (dn) 2 . (4) 
where here, and in the remainder of this paper, we use the mostly plus metric convention. 

As detailed in [TJ], only certain choices of F in ^ will lead to theories that have second 
order equations of motion; the Lovelock invariants and their boundary terms. The terms 
available are 

£2 = -V=g , (5) 

£ 3 = V=gK , (6) 

C 4 = -^g~R, (7) 
3 

£5 = ^\ f -9^-GB , (8) 



where 



/C GB = --K% + KKl v - -K A - 2(R flu - -Rg, v )K^ (9) 

is the Myers boundary term from the second order Lovelock invariant in the bulk [22], and 
£3 is the Gibbons-Hawking- York boundary term for the Einstein-Hilbert action in the bulk 
[231 M- 

In the following, we use the notation II for the matrix of partials = d^d^ir, and [IP] = 
Tr(n n ), e.g. [II] = Dtt, [n 2 ] = d^Ttd^TT, as well as [vr n ] = dix ■ H n ~ 2 ■ drr, e.g. [tt 2 ] = 
d^ird^n, [tt 3 ] = d^Tid p d u Tid v Ti. Indices are raised and lowered with the flat metric, and we 
use the mostly plus signature. Also, we define 



In terms of these combinations of the field ir and its derivatives, the terms above become 



£ 2 = - v /1 + (9tt)2, (11) 

A = -[n] + 7 2 ^ 3 ] > (12) 

£ 4 = - 7 ([nf-[n 2 ])-2 7 3 ([vr 4 ]-[n][vr 3 ]) , (13) 

c 5 = - 7 2 (P 3 + 2 [n 3 ] - 3 [n] [n 2 ]) - 7 4 (e [n] [vr 4 ] - e [vr 5 ] - 3 ([n] 2 - [n 2 ]) [vr 3 ]) , 

(14) 

where we have explicitly retained all total derivatives. 

In 3 + 1 dimensions, the above terms are the only ones possessing the symmetry Q and 
yielding second order equations of motion. The first term is the DBI action, which when 
expanded gives the standard kinetic term for the scalar. The second is the relativistic version 
of the cubic DGP 7r-lagrangian (up to a total derivative). These are the DBI generalizations 
of the galileons studied in [13J. The galileons are recovered by expanding in powers of the 
field 7r and taking the lowest non-trivial contribution from each term 



The resulting equations of motion take the form S n = 0, with n — 2, 3, 4, 5, and 

S 2 = 7 [n] - 7 3 [tt 3 ] , (15) 

£ 3 =7 2 ([n] 2 -[n 2 ])+27 4 ([vr 4 ]-[n][vr 3 ]) , (16) 

S A = 7 3 ([II] 3 + 2 [IT 3 ] - 3 [IT] [n 2 ]) + 7 5 (6 [IT] [vr 4 ] - 6 [vr 5 ] - 3 ([IT] 2 - [IT 2 ]) [tt 3 ]) , 

(17) 

s 5 = 7 e (jn] 4 - e [n 2 ] [n] 2 + 8 [n] [n 3 ] + 3 [n 2 ] 2 - 6 [n 4 ] ) . (is) 

These satisfy the following interesting recursion relation noticed in [14], 

JiW=9)-K t (19) 

±-(^K)=R, (20) 

^W=jR) = \KQ B , (21) 
6 2 

■^W=9Kgb) = 3>C GB4 , (22) 
where Cgb 4 = R 2 — 4-R 2 ^ + R 2 ^^ is the second order Lovelock invariant. 
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In this paper we consider a theory containing these terms with arbitrary coefficients d n , and 
which is linearly coupled to the trace T of the energy momentum tensor of matter, so that 
the complete Lagrangian density is 

5 

C = ^d n C n + nT , (23) 

n=2 

with equation of motion £ = 0, where 

5 

S = Y^d n S n + T. (24) 

n=2 

The linear coupling is not invariant under the symmetry operation (fl]). Rather, it was chosen 
for simplicity and for comparison with the results of [13] where the same choice was made. 
It is also the coupling that arises if the scalar is considered as a modification to gravity 
that conformally mixes with the graviton, as happens in the DGP model. Although there 
may exist physically interesting couplings which obey the symmetry, the simplest example, 
d fJ- 7rd u 7rT' J ' u ', which arises naturally from the brane construction, gives no contribution to the 
equations of motion for static sources. 

Our goal is to derive constraints on these models from the requirements of stability and 
subluminality of mode propagation around spherically symmetric backgrounds. We shall 
begin this analysis in the next section, but it is important to note that one constraint can 
be seen immediately; 

d 2 > , (25) 
since otherwise the kinetic term will yield a ghost (or will be absent, if we set d 2 = 0). 



III. SPHERICAL SOLUTIONS 



We search for static spherically symmetric solutions to the equations of motion in spherical 
polar coordinates (r,9,(f>), in the presence of a positive mass delta function source at the 
origin 

T = -M5 3 (r), M > . (26) 

To evaluate the equations of motion we need find only the non-vanishing elements of = 
d^d u TT — T^ u d a ii. These are Ii rr = ir irr , Ugg = r7r r , and = rsin 2 ^^^. Since the flat 
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metric is diagonal we then have 

[IP] = (n rr ^)"+(n^) n +(n^) n 



2% n r 

7T„„ H 



\n+l 2 n 

= 7T r 7T rr 



(27) 
(28) 



Using these, the equations of motion (24) become 



£2 = 


1 d r 






r 3 y 


£3 = 


2 d r 




^[ 


r 3 y 


£4 = 


2 d r 




^[ 


r 3 y 


£5 = 


. 





where we have defined 



y 



r 



(29) 
(30) 
(31) 
(32) 

(33) 



The fifth order term vanishes because our focus on static solutions reduces the problem to a 
three dimensional one, and the fifth order term is trivial in three dimensions. The remaining 
equations of motion can be written as a polynomial in y as 



1 d 

r 2 dr 



[r 3 P{y)] = M5 3 (r) 



with 



P(y) = d 2 y + 2d 3 y 2 + 2d 4 y 3 



(34) 



(35) 



Note that the equations of motion are a total r-derivative. This is a consequence of the 
shift invariance 71 — > 71 + c of the Lagrangian, which has an associated Noether current J M , 
in terms of which the equations of motion take the form c^(— J M ) = 0. We may therefore 
integrate the equations of motion once to obtain 

M 



P(y) 



Anr 3 



(36) 



We now study the existence of spherically symmetric solutions, and the resulting constraints 
on the coefficients ^2,^3,^4- Our boundary condition is that it approaches a constant as 
r — > 00. The other boundary condition is fixed by the delta function at the origin. Focusing 



on small r, (36) yields 



7T 



13 



75-/2)3/2 

7 



d.i 



M 
8^ 



(37) 



This determines a finite value for it' at the origin, and therefore implies that 7r must also be 
finite there. Since the absolute value of the prefactor in front of dj on the left hand side is 
always less than unity, we then obtain the constraint 

r , . M 

\ck\ > — . (38) 

This constraint is unique to the DBI action - no such constraint arises in the usual galileon 
theories. The fourth-order term dominates at short distances, and its non-linearities render 
7r finite at the origin. In particular therefore, note that there are no spherically symmetric 
static solutions in the pure DBI model, for which = = 0. 

As we have demonstrated, 7r'(r) ranges from some finite non-zero value at r = 0, to zero as 
r — > oo (since 7r itself goes to a constant). Thus, the variable y = jtt' jr ranges from infinity 
to zero as r ranges from zero to infinity (we will see shortly that it does so monotonically) . 



As r varies from the origin to infinity, the right hand side of ( 36 ) ranges from zero to infinity, 
so the cubic polynomial on the left must do so as well. Looking at small y, along with the 
requirement d 2 > for a healthy kinetic term, tells us that P(y) intersects the origin and is 
monotonically increasing near the origin, and hence that y as a function of r is monotonically 
decreasing in the same region. As y gets larger (r smaller, P(y) larger), the solution for 
y(r) must continue to exist and be smooth, which means that P(y) must not have any of 
its critical points in the region y > 0. Thus P(y) monotonically increases for y > 0, and 
hence y(r) is monotonically decreasing for r > 0. Looking at the form of y, this implies in 
turn that ir'(r) is monotonic, ranging from some finite value to zero as r goes from zero to 
infinity. Integrating, we see that 7r(r) is monotonic as well. 

The condition we have then is 

P'( y ) =d 2 + 4:d 3 y + 6d 4 y 2 > 0, for y > . (39) 
Focusing on large y implies that d 4 > 0, so that we can now remove the absolute value sign 



in (38). We already know that d 2 > 0, from the requirement of a healthy kinetic term, but 



it is worth pointing out that a direct implication of (39), applied at small y, is that spherical 
solutions do not exist for a ghost-like theory with d 2 < 0. Furthermore, we are safe if the 
minimum of P'(y) occurs above zero, which happens if 



\j 2 d-2<h ■ (40) 



Otherwise, the largest root of P'{y) must occur for y < 0, which happens if d 3 > 0. 

In summary, the flat space theory is ghost-free and spherical solutions exist if and only if 



d 2 > 0, d 4 > d 3 >-J^d 2 d 4 . (41) 



IV. STABILITY 

The existence of spherically symmetric solutions is, of course, not sufficient to guarantee 
viability of the theories in question. The next test is to examine the stability of these 
solutions. To do this, we expand the action in perturbations around the spherical solutions 

ir(x) = n (r) + ip(x) , (42) 

and isolate the terms quadratic in <p. These terms take the form 

S ^ = lJ dt J d2n J ^ (r)v " 2 ~ K r(r)(d r <p) 2 - K n (r)(d n <p) 2 ] , (43) 

where overdots denote time derivatives, (<9q<^) 2 = (deep) 2 + —^(^y?) 2 is the angular part of 
(Vy?) 2 , and the kinetic coefficients K depend on r through the background radial solution 
7r (r) and its derivatives. Note that the quadratic action contains only second derivatives 
acting on the perturbations. This is because the field equations are second order, despite 
the fact that the lagrangian is higher derivative, as we mentioned earlier. 

In order for the solution to be stable, each Ki(r) (i = t, r, Q) must be positive for all 
r > 0. If K t is negative in some region, then localized excitations will be ghostlike and will 
carry negative energy. If either of K r , Kn are negative in some region, then it is possible 
to find localized perturbations for which gradients lower the energy of the background solu- 
tion. This kind of instability, associated with negative gradient energy for certain classes of 
fluctuations, is more troublesome than a tachyon-like instability associated with a negative 
mass squared term or upside down potential. A tachyon-like instability is, like the Jeans 
instability, dominated by modes with momenta of order the tachyonic mass scale, which can 
be parametrically smaller than the UV cutoff, and thus computable within the effective the- 
ory. By contrast, the gradient instability can be due to very short wavelength wave-packets 
with high momentum. Thus, this instability also plagues fluctuations right down to the UV 



cutoff of the theory, so that quantities such as decay rates are dominated by the shortest 
distances in the theory, and cannot be reliably computed within the effective theory. 

To obtain explicit expressions for the functions Ki(r), we expand the equations of motion 
to linear order in (p 

£[ko + ^\^J^ = - R t( r ) <P + ^ d r{r 2 K r {r) d r p) + K a (r) d 2 Q p , (44) 
where <9q = ^g-§g ( sm #Jj?) + ^r^^i is the angular part of the laplacian. 
We begin with the radial perturbations, and find K r simply by perturbing the radial equa- 



tion (34), using a perturbation that depends only on r 



58 = i J- [r*P>(y)5y] = I J- [r 2 P\yWp'] . (45) 

From this we read off 

K r (r) = ^P\y) . (46) 



From (39), we then see that if the solution exists, then K(r) is automatically positive, since 
7 > 0. 



Now turn to the angular perturbations. To find Kq, we vary the full equations (24), al- 
lowing the perturbation to depend only on angular variables, and keeping in mind that the 
background depends only on r. Using the following useful expressions 



71TC 



/n-1 



5[U n }='-^ zr d 2 nV , 5[rc n ] = 0, <5 7 = 0, (47) 
it is simple to show that 

*n(r) = £± [r 2 P'(y)} . (48) 

Recall that the coefficient does not enter in either K r or Kq, because we are still consid- 
ering static configurations, for which the fifth DBI term vanishes. 

Lastly, we consider the temporal perturbations. We find K t by varying the full equa- 



tions (24), this time allowing the perturbation to depend only on time. Once again, some 



useful expressions 

5 [II] = -Cp , 5 [IT] = (n > 1) , 5 [n n ] = , 5 7 = , (49) 
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allow us to show that 



Kt(r) 



7 d 



[r 3 (d 2 + Qd 3 y + l8d 4 y 2 + 24d 5 y 3 )] . 



(50) 



3r 2 dr 

We see that d^ enters here for the first time, since we have deviated, at last, from static 
equations. 

dy 

As we have written them, the functions KAr) depend on 7, r, — and y. However, we 

dr 

dy 

may eliminate — in favor of y by using the implicit function theorem on the function 

dr 

M 

F(y, r) = P(y) — - _ „ = 0. This yields 



Airr 3 



dy d r F 3 P(y) 



dr d y F r P'(y) 
Substituting this into our expressions for the Ki(r) we obtain 

K r = 7 3 [da + 4d 3 2/ + Qd 4 y 2 } , 

~d\ + 2d 2 d 3 y + (Ad 2 3 - Qd 2 d A ) y 21 



K t = 1 



d 2 + 4:d 3 y + 6d 4 y 2 

d\ + {Ad 2 d 3 ) y + 12 (d 2 - d 2 d A ) y 2 + 24 (d 3 d 4 - 2d 5 d 2 ) y 3 + 12 (3d 2 - Ad 3 d b ) y A 

d 2 + 4d 3 y + 6d A y 2 



(51) 



Note that the explicit r dependence has canceled out. 



Since the solution spans all positive values of y as r varies from zero to infinity, we require 



K t and Kq, to be positive for all y > 0. The denominators in (51 ) are automatically positive, 



from (39). Given the constraints (41), The numerator in is positive for d 3 > J\d 2 d A 



which also ensures that the numerator in K t is positive provided d$ < 
The radial solution therefore exists and is stable if and only if 

d 3 > 



d 2 > , <i 4 > 



M 



in 



-dod. 



2«4 



3d 2 
Ad 3 



(52) 



V. PROPAGATION SPEED OF FLUCTUATIONS 



As a final test of the viability of the DBI galileon theories, we consider the propagation 
speeds of small fluctuations around the stable spherical solutions. For radially propagating 
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fluctuations this speed is 



K r 



At large distances from the source (small y), this becomes 

c 2 r = l + A + Otfl*) > 1 , 
where here and in what follows we express 7 in terms of y via 



(53) 



(54) 
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2/3 



2„2 



\4*P{y)J V 



(55) 



Therefore, given the constraints implied by existence and stability of the solutions, this is 
always super luminal. 



At smaller distances (larger y), the speed is 



3d 2 4 



3d\ — 4d 3 d 5 



M 
87rcL 



2/3" 



O 



(56) 



so the propagation speed is subluminal in this region if 

, 3dj ( M \ 2/3 
4 < 



44 V %nd/± 



(57) 



The speed of angular excitations is 



Ku 
K t 



(58) 



The difference between the numerator and the denominator is, apart from an overall positive 
factor, 

K n -K t ~ -2d 2 d 3 y - (8dl - Qd 2 d 4 ) y 2 - 24 {d 3 d 4 - 2d 2 d 5 ) y 3 - 12 (3d 2 - Ad 3 d 5 ) y 4 . (59) 



Given the constraints (52), this is always negative, so the speed of angular excitations is 



always subluminal. Also, the angular speed goes to zero as r goes to zero. The radial and 
angular speeds for a sample solution are shown in figure [TJ 

Certainly the existence of superluminally propagating modes raises questions about the via- 
bility of galileon DBI theories. Whether such a feature in really a problem that conclusively 
rules out a low-energy effective theory is still being debated [25H27] , but it has been argued 
that, at the least, it may preclude the possibility of embedding the theory into a local, 
Lorentz invariant UV completion [28J. 
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FIG. 1: Speed of fluctuations and Cq, in the radial and angular directions respectively, for a 
sample solution satisfying the existence and stability constraints (52), as well as (57). The values 
chosen are cfo = 1, d% = 2, = 1, d$ = —1, M = 1. 



VI. CONCLUSIONS 



The DBI galileon theories establish a natural generalization of, and connection between, the 
galiieon and DBI models through their higher-dimensional realizations and brane actions. 
In this paper we have studied spherically symmetric solutions to the DBI galileon models, 
demonstrating that there exists a range of parameters in which such solutions exist. We 
have also examined the stability of these solutions and computed the propagation speeds 
of perturbations around the solutions. While we have found that there exists a region 
of parameter space in which our solutions are stable, we have shown that these solutions 
always exhibit superluminal propagation. Such behavior is familiar from that of the ordinary 
galileon theories. Thus, although one might have thought that the 7 factors appearing for 
DBI galileons could cure the superluminality issues, the results we find here indicate that 
they do not. 

We have worked in dimensionless units, which corresponds to setting to unity a scale, A, 
suppressing all the non-linearities in the Lagrangian. In addition, we have absorbed into 
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the stress tensor a scale, M p , representing the coupling strength. Restoring these scales, the 



condition (38) tells us d 4 > M/M p , so in gravitational applications, where M is the mass of 
the Sun and M p the Planck mass, this tells us that <i 4 must be huge, of order the solar mass 
in Planck units. One might worry that this necessitates strong coupling, but this is not the 
case because the coefficient d%, which multiplies the kinetic term, may also be chosen to be 
very large, so that after canonical normalization the true couplings are still small. 

To see the consequences of this, consider expanding the action with the scale A restored. The 



DBI term reads schematically c^A^^l -|- ^4 ~ (d%) 2 + -^(dfc) 4 H , with the canoni- 

cally normalized field 7r = d^ 2 ir. The scale suppressing the non-linear terms here is d^A. 
Similarly, the quartic galileon term is, schematically, c? 4 1 + + • • • ^-(<9 2 7r) 2 (<97r) 2 = 

^4 (^ 2 7r) 2 (^7r) 2 + ^2 dhs (<9 2 7r) 2 (<97r) 4 H , which means that the strong coupling scales 

are [ 1 A, ( ^ J A, • • • . Since <i 4 is so large, keeping the lowest strong coupling scale 
reasonably high requires choosing g? 2 large, say d\ ~ <i 4 , in which case all the higher order 
DBI scales are much higher (corresponding to small coupling), and the theory becomes very 
similar to the ordinary galileons, explaining why we find conclusions similar to the con- 
clusions in that case. In addition, note that the coupling to the stress tensor, in terms of 
the canonically normalized field, is ~ -575 — 7tT. so that the true Planck mass is actually 

d 2 M p 

1/2 

~ d 2 M p , and the necessary size of d 4 is actually larger than the solar mass in physical 
Planck units. 

On the other hand, in some situations, it may be too much to demand that the spherical 
solutions exist for all r. For example, if n represents the fifth coordinate of a brane embed- 
ding, we should not expect that the brane configuration should be everywhere expressible 
as a single valued function of the four coordinates x M (the solutions of |2HrE2] are examples 
of this). In this case, the restrictions on the coefficient o? 4 may be relaxed. 

DBI galileon theories therefore, like the ordinary galileons, face a challenge from the superlu- 
minal propagation of perturbations around simple spherically symmetric solutions. Whether 
these theories are viable depends on the development of an argument that this superluminal- 
ity does not lead to the pathologies that are traditionally associated with this behavior, or 
whether a modification to the theory or its couplings to matter or gravity can eliminate this 
behavior. If the effects of gravity are taken into account, It should be mentioned that the 
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coupling of galileons to gravity is non-trivial if one wishes to keep the equations of motion 
second order [331 EH] , and the issue of superluminality should in principle be re-examined 
in the full covariant context, though the effects should be Planck suppressed. It should 
be straightforward to extend these results to anti-DBI type theories [351 EE] , which do not 
necessarily have a higher dimensional brane interpretation. 
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